Phase diagram of glassy systems in an external field 
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We study the mean-field phase diagram of glassy systems in a field pointing in the direction of a 
metastable state. We find competition among a "magnetized" and a "disordered" phase, that are 
separated by a coexistence line as in ordinary first order phase transitions. The coexistence line 
terminates in a critical point, which in principle can be observed in numerical simulations of glassy 
models. 



The crossover from supercooled liquids to glasses, 
although manifestly an off-equilibrium phenomenon, 
presents characteristics reminiscent of equilibrium phase 
transitions Probed at fixed observation time scale, 
the various thermodynamic quantities show non-analytic 
behavior at the glassy temperature Tg. While the singu- 
larity is pushed to lower temperatures for increasing ob- 
servation times, it does not seen to decrease in intensity. 
Whether for infinite probing times the singularity persists 
to produce a "true" finite temperature phase transition is 
a longly debated question, and is at present experimen- 
tally and theoretically unsolved. Experimentally, due to 
the rapid growth of the relaxation time as temperature is 
lowered, it is difficult to equilibrate the systems at tem- 
peratures that exceed the estimated critical temperature 
by less then 15 %. Theoretically, the only well estab- 
lished issues rely on mean-field approximation. Mean 
field theory for disordered models that are conjectured 
to be in the same universality class of structural glasses, 
predicts a thermodynamical transition, along the Gibbs- 
DiMarzio scenario of the entropy crisis at finite tempera- 
ture Tg . The transition at Tg combines features of first 
order and second order phase transitions. Like in a sec- 
ond order phase transition the energy is continuous and 
the specific heat has a jump. Like in a first order phase 
transition there is no divergent susceptibility. Detailed 
studies of the mentioned models Q have shown that it 
exists an interval of temperatures Tg < T < Td where 
the thermodynamic, as well as the dynamic is dominated 
by the existence metastable states, which can trap the 
system for infinite time. In finite dimension the picture 
must necessarily be modified, as metastable states can 
not have a infinite time life. An appealing possibility is 
that metastable states continue to exist with finite decay 
time. Transitions among metastable states would domi- 
nate the dynamics in that region. The inclusion of these 
phenomena is one of the most challenging problem in the 
theory of glasses [D . 

In ordinary mean-field theories the appearance of 
metastable states is signaled by local minima of the "po- 
tential function" , the free-energy as a function of the or- 



der parameter, and can be stabilized by the introduc- 
tion of an external field coupled with the order param- 
eter. In glasses, while the system is dynamically frozen, 
long range order is absent and it is impossible to iden- 
tify any intrinsic (physical) order parameter for a single 
system allowing to distinguish different thermodynami- 
cal states. The usual Edwards- Anderson parameter, that 
measure the degree of freezing, takes the same value in all 
the dominant states. Spin glass theory jj] suggests that 
measures of distance (or alternatively similarity) among 
different "real replicas" can be taken as appropriate or- 
der parameter. A potential, or better a class of potential 
functions, can be defined considering the partition func- 
tion of several real replicas with constraints on the mu- 
tual distances Q . Alternatively one can release the con- 
straints and introduce a field coupled to the distances. 
The two constructions are related by the Legendre trans- 
form, and both have been employed to study glassy sys- 
tems |-|. This point of view, although it can not help in 
real experiments, can give theoretical insights, and can be 
useful in numerical experiments. Different question can 
be addressed with different potential function. In this 
paper we limit ourselves to potentials involving two real 
replicas. More complicated three-replica potential have 
been considered in |^Jl^. Let us, to be defined, consider 
a system of N interacting variables Si, i = 1, N , with 
Hamiltonian H{S), and be q{S,S') an overlap function, 
i.e. an intensive measure of similarity among the con- 
figurations S and S'. We can consider two replicas of 
the system trying to thermalize at the same temperature 
with an attractive coupling. The appropriate thermo- 
dynamical description of this situation is given by the 
partition sum: 

Za^Y. exp[-/?(iJ(Si) + i7(S2))+/3e7Vg(S,S')]. 



Si ,82 



(1) 



We will call r^(e) — — T/N log Z a the annealed poten- 
tial. Alternatively, we can consider a reference config- 
uration a thermalized at temperature T', acting as an 
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external field on a replica S which tries to thermalize at 
temperature T. In this case the right partition function 
is 

ZQ=^cxp[-/3i/(S)+/5eiV9(S,a)], (2) 
s 

and rQ(e) = — T/A^logZg will be the quenched poten- 
tial. The idea of the quenched potential is particularly 
appealing for the physics of glasses. At a given cool- 
ing rate, the system equilibrates within the super-cooled 
liquid phase, until the glassy transition temperature Tg 
is reached. Below that temperature it is reasonable to 
think that the system will be stacked for a long time in 
a metastable state. While fast degree of freedom ther- 
malize at temperature T, the slow degree of freedom 
will remain confined for a long time close to the equi- 
librium configuration reached at Tg. In such conditions, 
restricted Boltzmann-Gibbs averages generated by a par- 
tition function of the kind with T' identified with Tg 
could be appropriate. In this pseudo-thermodynamic de- 
scription the parameter e, which fixes the equilibrium 
value of the overlap with the reference configuration, 
should be considered as a dynamical variable, fixed by the 
condition that the system has equilibrated in the neigh- 
bors of the reference configuration cr, but has not had 
the time necessary to escape from the metastable state. 
Within the limit of validity of this picture the usual laws 
of thermodynamics must hold in the glassy phase. The 
picture certainly loose its validity for times large enough 
to escape from the metastable state and aging dynamics 
sets in the system 

The quenched case in an infinitesimal field and T — T' 
was studied in detail in [Q. In this paper we extend 
the analysis to the equilibrium phase diagram of the 
quenched potential in the T — e plane. For definite- 
ness let us introduce the model that we will use in 
the following, the same lines of reasoning can be fol- 
lowed in general. The model we study is a spherical 
p-spin model [^. This is defined in terms of N real 
dynamical variables St, [i — l,...,iV) subjected to the 
constraint X^iLi — ^ ^'^'^ interacting via a random 
Gaussian Hamiltonian H{S), defined by its correlator 
H{S) H{S') = Nfiqs^s')- The overlap qs,s' is a mea- 
sure of similarity among configurations and is given by 
qs,s' = ^I^Tlii^i^i- If the correlation function is taken 
of the monomial form J{q) = l/2qP the Hamiltonian 
can be represented as H = -J2]l<...<ip Jii,...,ipSi^...Si^ 
with independent centered Gaussian couplings Jii^...^i^, 
which motivates the name of the model. We will stick to 
this case throughout all this paper. The model has been 
studied extensively during the last years, and furnishes, 
for p > 2, a good mean-field model of fragile glass. It 
has been often observed that the Langevin relaxation 
of this model leads to equations homologous to those of 
mode coupling theory [Q. 



The quenched free energy TQ[a\[3,e\ should be self- 
averaging with 
respect to the distribution of cr, so we can compute it 
as r(e,T,r') = H{ct))T /N\ogZ[a-e]. 

In the following we will study the phase diagram of the 
model in the e— T plane in two situations: a) T' ~ T, cor- 
responding to restricting the partition sum to the vicinity 
of a particular equilibrium state at each temperature, b) 
T' fixed, corresponding to probe the evolution of the free- 
energy landscape in the vicinity of a fixed configuration 
of equilibrium at T' when T is changed. 

The Legendre transform of r(e,T,r'), V{q,T,T') = 
miug T{e) + eq, which corresponds physically to constrain- 
ing the value of gs.o- to q, was studied in detail in |^ with 
the aid of the replica method. The interested reader can 
found there details about the general method and the an- 
alytic expression for V in the case of the spherical p-spin 
model. 
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FIG. 1. The potential as a function of q for 

T' = 2TsTd/{Ts + Td) and various values of T, in order of 
decreasing temperatures from top to bottom. 

The shape of the function V turned out to be the 
characteristic one of a mean-field system undergoing a 
first order phase transition. At high enough tempera- 
ture V is an increasing and convex function of q with a 
single minimum for q — Q. Decreasing the temperature 
to a value T/, where for the first time a point qf with 
^"(9/) = appears, the potential looses the convex- 
ity property and a phase transition can be induced by a 
field. A secondary minimum develops at Td, the temper- 
ature of dynamical transition |^ , signaling the presence 
of long-life metastable states. The minimum of the po- 
tential has received a dynamical interpretation as 
corresponding to the states reached at long times by the 
evolution at temperature T starting at time zero from 
an equilibrium configuration at temperature T' . This al- 
lows to fix q in the pseudo-thermodynamical description 
mentioned above to the value it takes in the minimum 
of the potential. The height of the secondary minimum 
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reaches the one of the primary minimum at T = Tg and 
coexistence in zero field takes place. This is the usual 
statical transition point in zero field, and it is not ac- 
companied by the release of latent heat. In figure 1 we 
show the shape of the potential in the various regions. 
The attentive reader would have noticed that respect to 
the curves presented in only one secondary minimum 
is present at low temperature. The results we present 
here are corrected taking into account replica symmetry 
breaking effects, the meaning of which will be discussed 
in . Although the behavior of the potential function 
is analogous to the one found in ordinary systems under- 
going a first order phase transition the interpretation is 
here radically different. While in ordinary cases different 
minima represent qualitatively different thermodynami- 
cal states (e.g. gas and liquid), this is not the case in the 
potential discussed here. In our problem the local mini- 
mum appears when ergodicity is broken, and the config- 
uration space splits into an exponentially large number 
of components. The two minima are different manifesta- 
tions of states with the same characteristics. The height 
of the secondary minimum, relative to the one at g = 
measures the free-energy loss to keep the system in the 
same component of the quenched one. At equal temper- 
atures T = T' this is just the complexity TS, i.e. the 
logarithm of the number of distinct pure states of the 
system. For T 7^ T' it also takes into account the free- 
energy variation of the equilibrium state at temperature 
T' when followed at temperature T. 
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FIG. 2. Phase diagram in the e - T plane for T' = T. 
The upper curve is the spinodal of the disordered state, the 
lower one the spinodal of the magnetized state, and the middle 
curve the coexistence line. The coexistence line touches the 
axes e = at r = Tji, while the spinodal of the magnetized 
state touches it at T = T — d. For T <Ts the spinodal of the 
disordered state remains finite and touches the T — axes at 
finite e. 

The presence of the field e adds finite stability to the 
metastable state, and the transition is displaced at higher 
temperatures. In figure 2 we display the phase diagram 



of the p = 4 model in the case T' = T. The coexistence 
line departs from the axes e = at the transition tem- 
perature Ts and reaches monotonically a critical point 
{Tcr, £cr)- In figure 2 it is also shown the spinodal of the 
"magnetized" solution, which touches the e = axes at 
the dynamical temperature Td, and the spinodal of the 
"disordered" solution for temperaratures larger then T^. 

The coexistence line for T' fixed in the interval Td < 
T' < Ts is qualitatively similar to the one of figure 2 
at high enough temperature, but (for T' > T^) it never 
touches the axes e = 0. 
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FIG. 3. Phase diagram in the e — T plane for 

T' = 2TsTa/{Ts + Ta). The upper curve is the spinodal of 
the disordered state, the lower one the spinodal of the mag- 
netized state, and the middle curve the coexistence line. 

Even at zero temperature there is a first order phase 
transition in e, refiecting the fact that the ground state of 
the system is lower then the energy of the reference state 
when followed at T = 0. This can be seen in figure 3 
where we show the phase diagram for T' = 2TsTd/ {Ts + 
Td). At the critical point the transition is second order. 
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FIG. 4. Latent heat of the transition as a function of the 
temperature for T' as in fig. 3. 
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While the transition in zero field is not accompanied 
by heat release, a latent heat is present in non zero e. 
In figure 4 we show, in the same conditions of fig. 3, the 
latent heat Q = — — e{q+ — q^) where E+ {q+) and 
E- are the internal energies (overlaps) respectively 
of the magnetized and unmagnetized states. Notice that 
(as it should) the latent heat is zero at the critical point 
and at T = 0. The magnetized state roughly reflects 
the properties of the equilibrium states at temperature 
T' followed at temperature T, while the disordered state 
reflects the properties of the true equilibrium states at 
temperature T. We see that at high temperature the 
magnetized state is energetically favored, while at low 
temperature it has an energy higher then the one of equi- 
librium. The point where Q changes sign does not cor- 
respond to a second order phase transition. Finally, in 
figure 5 we show, for a fixed temperature the curve of 
q{e) obtained by the Maxwell construction. 
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FIG. 5. Equation of state for T' = 2TsTd/{Ts + Ta) and 
T — 0.609. The horizontal line corresponds to coexistence 
and is obtained by the Maxwell construction. 

So far we have presented results for the quenched po- 
tential Fq. In the annealed case, it was shown in Q the 
potential Va has qualitatively the same shape than the 
quenched one (Vq); it can therefore be expected a phase 
diagram similar to the one of figure 2 also in this case. 

Although we have based our discussion on a mean- 
field model, we expect that the qualitative features of 
the phase diagrams presented survive in finite dimension. 
We belove that the existence of a coexistence line, termi- 
nating in a critical point, is a constitutive feature of sys- 
tems whose physics is dominated by the existence of long 
lived metastable states. The predictions of this paper can 
be submitted to numerical test in glassy model systems 
as like e.g. Lennard- Jones or hard spheres, or polymer 
glasses. For example the identification of the complexity 
S as the free energy difference between the stable and the 
metastable phases allows a direct measure of this quan- 
tity in a simulation. Indeed the ending of the transition 



lines in a critical point implies that the metastable state 
can be reached via closed paths in phase diagram leaving 
always the system in (stable or metastable) equilibrium; 
and the free energy difference of the two phases computed 
integrating the specific heat a long the loop. 
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